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Recently, it has been shown that the quantum equilibrium distribution in the original Bohm’s 
model is unstable and so it isn’t a tenable physical theory [Proc. R. Soc. A 470 20140288 (2014)]. 
In this paper we show that a natural modification of the Bohm’s quantum force leads to a stable 
quantum equilibrium without any change to statistical predictions of this model in equilibrium state. 
Moreover, it is shown that an initial non-equilibrium phase space distribution relaxes to equilibrium 
distribution via a coarse grained H-theorem. But before end of the relaxation process, for example 
in the early universe, this modihed quantum force can lead to the new predictions beyond standard 
quantum mechanics. 
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1. Introduction-. Bohmian mechanics is a completely 
coherent casual theory for describing quantum systems 
which introduced by L. de Broglie in 1927 [T] and re¬ 
formulated and developed by D. Bohm in 1952 [ 2 ]. In 
this theory, unlike the orthodox quantum mechanics, 
the quantum particles have an actual position in each 
moment and moving along well defined continuous tra¬ 
jectories under effect of the ^/’-field. In non-relativistic 
Bohmian mechanics, the time evolution of the wave func¬ 
tion, for a system of N spineless particles, is given by 
Schrodinger equation: 
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Where rrii is the mass of i-th particle and x = 
[xi, ...,xm) G But for describing particles’ motion, 

there are two approaches [SE]: de Broglie’s first-order 
dynamics and Bohm’s second-order dynamics. In the de 
Broglie’s dynamics, the state of the system is completely 
determined by the wave function ip and particles’ posi¬ 
tions X — {Xi ,..., Xn), while the law of particles motion 
is given by ” guiding equation”: 


( 2 ) 

Where S is the phase of the wave function, tp = 
and m is defined as m = diag(mi,..., mAr) and V = 
(Vi,..., Vjv). But in the Bohm’s reformulation of the 
de Broglie’s model, the state of the system is determined 
by wave function, particles’ position and also particles’ 
velocity. In this model, the particles’ law of motion is 
given by Newton’s law with an additional ’’quantum po¬ 
tential” : 

^^^ = -^iV + Q)L=xitv ( 3 ) 
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Where Q is defined as Q = (—^^/2mi)Vi^i?/i?. 

Both of the de Broglie’s and Bohm’s models are deter¬ 
ministic theories and so the ’’probability concept” is not 
considered as an intrinsic property in these models: i.e. 
the stochastic characteristics of the quantum phenom¬ 
ena must be explained by the aid of the stochastic initial 
conditions mm- It is well known that the empirical pre¬ 
dictions of the standard quantum mechanics follow from 
the de Broglie’s dynamics, if it is assumed that the initial 
positions of particles is distributed as 

p{x-,to) = \tp\^. (4) 

The distribution p = \ip\'^ is called ’’quantum equilib¬ 
rium distribution”. It is an elementary consequence of 
the de Broglie’s dynamics that the quantum equilibrium 
distribution is preserved in time: if p = \%p\'^ at to then 
p = \ip\'^ at all future times. This property is referred to 
as ” equivariance” [5]. Similarly, if the initial distribution 
of particles in phase space be as 

fBix,p;to) = R‘^S{p-VS), (5) 

Where p = (pi, ...,p 7 v), then the empirical predictions of 
the standard quantum mechanics follow from the Bohm’s 
dynamics. We name fB{x,p) = R‘^6{p — XS) as ’’Bohm’s 
quantum equilibrium distribution”. The Bohm’s quan¬ 
tum equilibrium is preserved in time in the Bohm’s dy¬ 
namics same as p = \'ip\'^ in the de Broglie’s dynamics 
(equivariance) [5]. The equivariance of the Bohm’s quan¬ 
tum equilibrium /b (x, p), ensures the validity of the Born 
rule and guiding equation in all times. So in the quan¬ 
tum equilibrium conditions, both of the models lead to 
the same trajectories and so are equivalent. But it should 
be noted that in the both formulations of the Bohmian 
mechanics, the wave function ip is considered as an actual 
field in the configuration space which affects on the parti¬ 
cle’s motion, so at the fundamental level, ip is not a priori 
related to probability notion, even through initial condi¬ 
tion. Firstly, W. E. Pauli objected that in a fundamen¬ 
tally deterministic theory, taking a particular distribu¬ 
tion p = \ip\‘^ as an initial condition is not reasonable [5]. 
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Notwithstanding this objection, D. Bohm had already 
suggested in his original papers that the distribution 
p = I'f/jp was similar to thermal equilibrium distribution 
in ordinary statistical mechanics, and must be derived by 
suitable statistical-mechanical arguments I31S]. In this 
regard, in 1991, A.Valentini showed in the de Broglie’s 
dynamics, an arbitrary non-equilibrium distributions re¬ 
lax to quantum equilibrium distribution, p —> based 

on a ’’sub-quantum H-theorem” Moreover, in re¬ 

cent years, the numerical simulations have shown that 
an initial non-equilibrium distribution p ^ rapidly 
approaches to \’tp\^ on a coarse-grained level [TMT5] . But 
the concept of ” non-equilibrium state” in the Bohm’s dy¬ 
namics is wider than the de Broglie dynamics; because 
in the former, the particles’ velocities may be deviate 
from VS'/to, beside the fact that p may be deviate from 
The question that arises at this point is whether or 
not the non-equilibrium distributions in Bohm’s dynam¬ 
ics, / ^ R^5{p— VS”), tend to relax towards equilibrium 
distribution fs- Because of two following facts, an ini¬ 
tial distribution with non-zero Lebesgue measure support 
will not be able to relax to Bohm’s quantum equilibrium 
[14| : I) The support of the Bohm’s quantum equilibrium, 
has zero Lebesgue measure in phase space, II) The phase- 
space volume is conserved in the Bohm’s dynamics. In 
addition, recently it has been shown that the quantum 
equilibrium in the Bohm’s dynamics - unlike de Broglie 
dynamics -is unstable i.e. a small perturbation to the 
equilibrium rapidly grows with time |I4j . To emphasize 
this instability, S. Goldstein and W. Struyve, showed that 
for various wave functions, the particles tend to escape 
the wave packet in Bohm’s dynamics m- Based on this 
fact, it has been concluded that the Bohm’s dynamics 
is untenable as a physical theory M- But it should be 
noted that the D. Bohm did not consider the original ver¬ 
sion of his model as a final theory and even suggested a 
general form for modifying the quantum force as follows 

m- 

m^^^ = -W{V + Q)+F{p-WS), ( 6 ) 

however, he did not introduce any explicit form for the 
function F {p — VS). In this paper, we derive a natural 
generalized quantum force that leads to a stable casual 
quantum theory. Moreover, we will show that in this 
modified Bohm’s dynamics, a non-equilibrium distribu¬ 
tion relaxes to equilibrium distribution as same as de 
Broglie’s first-order dynamics. 

2. Modification of the Bohm’s dynamics-. Since the 
mathematical origin of the instability of the Bohm’s equi¬ 
librium distribution is rooted in its singularity, we suggest 
replacing the Dirac delta function in eq.(5) by a Gaussian 
distribution with a width of p > 0 around VS : 

ffj. {x, p; t) = A^exp ^. (7) 

-h 

Notice that, the above ’’modified equilibrium distribu¬ 
tion”, /^, as same as the Bohm’s equilibrium distribution 


fs, is consistent with standard probabilistic interpreta¬ 
tion of the wave function, because: 

P = y//.(a;,p)d^^p= |^/’|^ (8) 

>/ = / {p/m)f^ {x,p) d^^p = \iljfVS/m. (9) 

It is clear that, is not equivariant under Bohm’s dy¬ 
namics, so if we want to consider it as a quantum equi¬ 
librium distribution then the Bohm’s dynamics must be 
modified. For this purpose, we derive the time evolution 
of the /^, using Schrodinger equation: 

Where Vp = (Vp^,..., Vp^^) is gradient operator in the 
momentum space and Fq is: 

Fq = -VQ-^^ + -VV{IS){VS-p), (II) 
mum 

where I is identity matrix in 3D position space. Equation 
(10) can be considered as a Liouville equation with an 
extra modifed quantum force Fq. So, we consider the 
particles’ equation of motion as: 

= {-VV + FQ)U^xit), ( 12 ) 

™ P=P(t) 

Where P = mX is particles momentum vector. This 
” modified Bohm’s dynamics” leads to the same Liouville 
equation as (10) for paticles distribution in phase space 
f{x,p;t). Therefor, if f=fp at some initial time tg then 
/=/p for all future times (equivariance). The equivari- 
ance of fp in phase space ensures the equivariance of I'i/’P 
in configuration space and based on the causal theory of 
measurement [4], the consistency between statistical pre¬ 
dictions of our model and standard quantum mechanics 
for all observables is guaranteed. However, in our model 
the average velocity of particles in each point is equal to 
VS/m , but the particles’ trajectories are completely dif¬ 
ferent from standard Bohemian trajectories, even if the 
initial velocities are chosen as X{tQ)=VS{X{tQ),tQ)/m. 
In fact, in equilibrium state, the standard Bohmian me¬ 
chanics could be considered as p, ^ 0 limit of our model: 
because in this limit, the modified quantum equilibrium 
distribution /p, reduces to the Bohm’s quantum equilib¬ 
rium distribution fs, and so we have P ^ VS for all 
particles in ensemble; Therefore in this case the modified 
quantum force reduces to standard Bohmian quantum 
force, Fq —VQ. 

As an example, consider the coherent states of a har¬ 
monic oscillator in one dimension, withF(a;) = kx^/2: 

= Nexp —5^((a: — acost)^ + i{t -\- 2xasint -\- a^sin2t)) 

( 13 ) 
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where N is normalization constant and a is an arbitrary 
constant. We choose h = k = m = 1, in this case, the 
eq.(12) leads to: 

X (t) = Xq ' ^ +cos(y7^t) (Xq — a) + a cos t (14) 

V k 

Where Xq and Xq are initial position and velocity of par¬ 
ticle respectively. In the limit /x —)• 0, the equation (14) 
reduces to corresponding equation of motion in Bohm’s 
dynamics: X {t) = Xq + XqI + a (cos t — 1). In this case, 
as in [20] has been mentioned, if Xq ^ l/m{dS{x, 0)/dx) 
Then the particle will escape to infinity and so the quan¬ 
tum equilibrium in Bohm’s dynamics is unstable. But 
even for a small amount of /i, our model circumvents this 
problem. In figure.I, particle trajectories in our model 
are compared with those of the Bohm’s model. It is 
clear that this result in not limited to the above example 
and our model generally avoids the problem; because the 
modified equilibrium distribution includes the parti¬ 
cles with p 7 ^ VS'. 

3. Sub-quantum H-theorem-. In this section, we 
want to show that the non-equilibrium distributions ap¬ 
proach to modified equilibrium distribution via an coarse- 
graining H-theorem similar Valentainis H-theorem. For 
this purpose, since our model is a second order dynamics, 
we must generalize Valentaini’s H-function, 

p In (p/IV’n , (15) 

from configuration space to phase space. For this, we 
introduce the counterpart quantities in the phase space 



FIG. 1. The particle trajectories in our model (with p = 1) is 
compared by Bohm’s trajectories in the case of the coherent 
state of harmonic oscillator (13), with a = 1. Wave packet is 
represented by the White shaded area. The Bohm’s trajecto¬ 
ries are represented by dashed curves and particles trajecto¬ 
ries in our model are represented by solid curves. All trajecto¬ 
ries have initial velocities as Xq = ±0.25 yf l/m{dS{x, Q)/dx) 
and the trajectories with same color have same initial condi¬ 
tions. The Bohm’s trajectories escape the wave packet but 
the trajectories of our model follow the wave packet. 


as: 


p{x) f{x,p), 

(16) 

IV’(a;)P ^ ft.ix,p). 

(17) 


So we redefine Valentinis H-function in our model as : 

H^ = Jdnf InC, (18) 

Where dft = d?^ x d^^p is differential volume element 
in the phase space and ^ is defined as i=f /This 
H-function is a non-negative quantity and will be zero 
if and only if f=ffj, everywhere, these features make it 
a useful measure of proximity to quantum equilibrium 
distribution. On the other hand, the ^ is preserved along 
particles trajectories in the phase space: 

|^(X(t),P(t);t)=0, (19) 

becuse / and obey identical Liouville equations in our 
model. In addition, fdfl is the number of the systems in 
ensemble which occupy a comoving volume dtd and so it is 
preserved along trajectories too. Therefore the exact 
is constant in time, dH^/dt=0, and so a non-equilibrium 
distribution can never really relaxes to equilibrium dis¬ 
tribution. But, since all physical measurements have a 
finite accuracy, it is reasonable to understand the re¬ 
laxation process in terms of a coarse-grained H-function 
which isn’t a conserved quantity, as same as classical case 
[Mini- For this purpose, we divide phase space into non 
overlapping cells of volume and define coarse-grained 
quantities as: 

and 

where / and are considered as constants in each cell. 
These values can be used to definitoin a course grained 
i7-function: 

H^= J dnj In?, (22) 

Where ? is defined as ?=///,i. Now, to prove that Hf^ 
decreases with time, same as the classical case, we assume 
no ’’microstructure” for the initial state, on the other 
word, we assume f{to)=f{to) and U{to)=ftJ.{to)- The 
proof now precedes along the lines of the Valentaini’s 
proof for first order Bohemian mechanics: From above 
assumption and the fact that the exact is constant in 
time, we have; 

(io) = {t) = I 


dn ff, ? In?. 


(23) 
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In addition, by using the fact that ^ is constant in each 
cell 6il, one can easily show that 

J dn , (24) 

and 

I an (? - c) = 0. (25) 


order dynamics and infact the only necessary properties 
for the modihed quantum equilibrium distribution are as 
follows: (I) It must be a local function of ip that satisfies 
equations (8), (9) and (II) Lebesgue measure of its sup¬ 
port sould be non-zero in phase space. As an example in 
definition of /^, one can use the Lorentzian representa¬ 
tion of the Dirac delta function instead of the Gaussian 
representation: 


So by using equations (23),(24) and (25), we have: 

= j dn /^(Cln(e/?) +?-e). (26) 


Finally, Since x ln(a;/y) -|- y — a; > 0 for all x, y, we have 
the subquantum H-theorem : 


dH^ 

dt 


< 0 . 


(27) 


The decrease of shows relaxation of non-equilibrium 
distribution to modified equilibrium distribution (/ —>■ 
Jfj). This coarse-grained H-theorem could be considered 
as the mathematical formulation of the Gibbs-Ehrenfest’s 
idea in the context of the second-order Bohmian mechan¬ 
ics: a non-equilibrium phase space distribution will tend 
to increase fine-grained microstructure and become more 
like equilibrium distribution on a coarse-grained level. It 
can be considered as a statistical derivation of the Born 
rule and also ” averaged guiding equation” in the context 
of the deterministic Bohm’s approach, without postulat¬ 
ing ” random fluctuations”, as done by Bohm and Vigier 
in 0. 

4 . Discussion and concluding remarks-. In this paper 
we introduce a modified quantum force which leads to a 
stable quantum equilibrium. But let us notice that the 
special form of the modified quantum equilibrium distri¬ 
bution /^, that we used in this paper is not the only pos¬ 
sible choice that could be used to make a stable second 




TT {p — V5')2 -I- ’ 


(28) 


which leads to a different form for quantum force via Li- 
ouville equation. In fact, there are many stable second 
order dynamics that lead to the same statistical predic¬ 
tions as standard quantum mechanics in the equilibrium 
state. But in principele for the non-equilibrium state, 
these modified dynamics may lead to different predictions 
that could be used as experimental test of them. Since 
all ordinary physical systems, such as atoms in the labo¬ 
ratory, have had a long astrophysical history, we can not 
see these systems in non-equilibrium state today. How¬ 
ever, before ending relaxation process in the early uni¬ 
verse, the non-equelibrium effects could be observed. In 
fact, it has been proposed that the GMB data may be 
used to set bounds on non-equilibrium deviations from 
quantum theory [I8h21j . As an out look, the CMB data 
can be used to determine some limits on the value of the 
parameter p. 
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